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ABSTRACT 

We investigate the luminosity dependent clustering of rest-frame UV selected galaxies at z ~ 4, 3, 
2.2, and 1.7 in the Keck Deep Fields (KDFs), which are complete to 1Z = 27 and cover 169 arcmin 2 . 
We find that at z ~ 4 and 3, UV-bright galaxies cluster more strongly than UV- faint ones, but at 
z ~ 2.2 and 1.7, the UV-bright galaxies are no longer the most strongly clustered. We derive mass 
estimates for objects in our sample by comparing our measurements to the predicted clustering of dark 
matter haloes in the Millennium Simulation. From these estimates, we infer relationships between 
halo mass and star formation rate (SFR), and find that the most massive dark matter haloes in our 
sample host galaxies with high SFRs (A/1700 < —20, or > 50 M© yr -1 ) at z ~ 3 and 4, moderate 
SFRs (-20 < Mi7 00 < -19, or - 20 M© yr" 1 ) at z ~ 2.2, and lower SFRs (-19 < Mi 700 < -18, 
or ~ 2 M© yr" 1 ) at z ~ 1.7. We believe our measurements may provide a new line of evidence for 
galaxy downsizing by extending that concept from stellar to halo mass. We also find that the objects 
with blue UV colors in our sample are much more strongly clustered than those with red UV colors, 
and we propose that this may be due to the presence of the 2175 A dust absorption bump in more 
massive halos, which contain the older stellar populations and dust needed to produce the feature. 
The relatively small area covered by the survey means that the absolute values of the correlation 
lengths and halo masses we derive are heavily dependent on the "integral constraint" correction, but 
the uniformly deep coverage across a large redshift interval allows us to detect several important 
trends that are independent of this correction. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
halos — large-scale structure of the Universe 



1. INTRODUCTION 

The study of high-redshift galaxies is a rapidly devel- 
oping field that is being driven in part by the ability to 
efficiently select large samples of galaxies using multi- 
wavelength photometry. The Lyman- break technique pi- 
oneered by Steidel and collaborators (jSteidel et al.l fl999. 
120031 [2004) working in the rest-frame UV has proven 
to be very successful at selecting star-forming galax- 
ies at z > 1.5, and at rest-fr ame optical wavele ngths, 
the BzK selection technique (|Daddi et al.l I2004T ) effec- 
tively locates galaxies with large stellar populations at 
1.4 < z < 2.5. Other high-z gala xy populations in- 
clude distant red galax ies (DRGs, iFranx et al.l 120031 : 
Ivan Dokkum et al1l2p03f). with J — K > 2. 3, extremely 
red objects (EROsT TThompson et al.l 119991) with large 
optical to near-i nfrared colors (R — K > 5), sub-mm 
galaxies (SMGs, iSmail et"aT1 I2004T). infrared- lum inous 
dust-obscured galaxies (DOGs HDev et ail [2008) with 
-P24 / Fr. > 100 0, and strong Lyman-a emitters (LAEs; 
IHu et al1 ll998). The large differences in selection crite- 
ria ensure that each sample represents a biased view of 
the total galaxy population at each redshift. Establish- 
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ing the properties of these different samples of galaxies, 
the relationships between them, and their evolution with 
redshift, is essential to our understanding of how galaxies 
form and evolve. 

On the theoretical side, it is well established that the 
clustering strength of dark matter haloes is strongly de- 
pendent on mass, with massive haloe s being the most 
strongly clustered (|Mo fc White! Il996t ). In our current 
understanding of galaxy evolution, the spatial distribu- 
tion of different types of galaxies is thought to represent 
differences in the mass distributions of the dark matter 
haloes that host them. For many samples of UV-selected 
high-z galaxies, an important link has been made be- 
tween clustering strength and UV magnitude (which is an 
indicator of star formation), wit h more luminous galax- 
ies clustering more strongly Ce.g-.lAdelberger et al.ll2005l : 
iLee et al.|[200l iKashikawa et al.l 12006! ) . This implies a 
direct link between halo mass and SFR at high-z, and 
it is important to refine this relationship in order to test 
models of star formation in young galaxies. Since there is 
no strong relation between SFR a nd clustering stren gth 
(halo mass) for galaxies at z ~ 0.3 (Hcinis et al. 20071) . it 
is important to identify the redshift interval over which 
a relation is in place, and also to determine if it extends 
to fainter galaxies. 

Clustering studies of high-z galaxies are no longer in 
their infancy, but many discrepancies and unanswered 
questions remain. One of the most significant current 
limitations in the study of UV-selected galaxy cluster- 
ing is the lack of uniformly deep UV coverage that 
probes significantly below the brightest objects at a large 
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range of redshifts. In particular, most large surveys, 
such as t he Great Observatorie s Origins Deep Survey 
f GOODS ■ iGlavarisco et al1l200l . lack deep U band cov- 
erage. In this respect, th e Keck Deep Fields (KDF, 
ISawicki fe Thompson! [2005) provide an important new 
perspective on galaxy clustering by allowing us to select 
galaxies with a broad range of UV luminosities (TZ = 23 
- 27) in four redshift samples at z ~ 4, 3, 2.2, and 1.7. 
The deep UV coverage of the KDFs over a large range of 
redshifts enables us to search for important trends that 
have not been previously reported. 

As in all of the papers in the KDF serie s, in the pr esent 
paper we use the AB flux normalization (Oke [1974ft and, 
unless otherwise noted, adopt fiM=0.3, J7a=0.7, Hq=70 
km s Mpc . To facilitate comparisons to previous 
studies, correlation lengths are expressed in units of h~ x 
Mpc, where h = 0.7. 

2. DATA 

In this section we briefly summarize some of the most 
important features of the KDF data on which we base our 
analysis. A detailed description of the KDF observations, 
d ata reductions, and high-z gal axy selection can be found 
in ISawicki fc Thompson! ([20051 ). while our z ~4 , 3, and 
2 luminosity functio n measurements are p r esented and 
discussed in detail in ISawicki fe Thompson! (|2006a|) . 

2.1. The KDF Survey 

Our analysis is based on results from our very deep 
multicolor Keck Deep Fields imaging survey carried out 
using a total of 71 hr of integration on the Keck I tele- 
scope. These KDFs use the same U n GlZI filter set and 
color-color selection techniques as a re used for b r ighter 
UV-selecte d galaxies in the w ork of [Steidel et al.l (I1999L 
20031 12004ft . In contrast to the lSteidel et all (jl999U2003l 
2004ft work, however, the KDFs reach TZu m =27 ; this is 
1.5 mag deeper than the Steidel et al. (1999, 2003, 2004) 
surveys and significantly below L* at z—2-4: even at z ~ 
4 we reach 2 mag fainter than M*. The KDFs have an 
area of 169 arcmin 2 divided into three spatially indepen- 
dent patches on the sky, with two ~7xl0 arcmin 2 fields 
and one ~7x5 arcmin 2 field. 

To Ku m =27, the KDFs contain 427, 1481, 2417, and 
2043 [/„G7?./-selected star-forming galaxie s at z ~4, 3, 
2.2, and 1.7, respe ctively, selected using the lSteidel et al.l 
(fl99l [20031 [2001 LBG, BX, and BM color-color selec- 
tion criteria. At our completeness limit, the KDF data 
probe galaxies with UV luminosities that correspond to 
SFRs (uncorrected for dust obscuration) of only 1-2 M Q 
yr _1 . Spectral energy distribution fitting suggests that 
there i s very little du st in galaxies at the faint end of this 
range (jSawickil 120111 ). so these star formation rates are 
low indeed. 

An important feature of the KDFs is that our use of 
the U n GlZI filter set lets us select high-z galaxies using 
the color-color selection criteria defi ned and spect r oscop- 
ically teste d to (TZ Hm = 25.5) by iSteidel et all (|1999l 
120031 12004ft . Our KDF sample is significantly deeper 
(Kiim = 27), making spectroscopic verification challeng- 
ing, but the commonality of filters makes us confident 
our selection is reliable (for a more detailed discussion, 
see ISawicki fe Thom pson 2005|). Nevertheless, we used a 
combination of published spectroscopy and new deep ob- 
servations to verify our selection to the extent possible. 



The KDF 09A and 09B fields ha ve been chosen to ov er- 
lap with the Q0933+2854 field of ISteidel et all (120031) in 
which these authors have spectroscopically observed a 
substantial number of their [7 n G7?.-selected LBGs, and 
secured redshifts for many of them. This overlap makes 
it possible to verify our selection, at least at the bright 
end of the KDF sample (TZ < 25.5). To this end, we 
cro ss-matched 19.0 < TZ < 25.5 KDF objects against 
the lSteidel et al.l (|2003ft catalog of C/ n G£-selected LBGs. 
The KDF 09A and 09B fields are almost entirely within 
the Steidel et al. Q0933+2854 field, with only 5% not 
overlapping (which we excluded from this analysis). In 
cross-matching, 120 of 139 Steidel objects are recovered 
in our survey. Several factors are likely responsible for 
us missing the remaining Steidel et al. LBGs: differ- 
ences in detection masks (e.g., to avoid pixels affected 
by bleeding and bright stars); photometric scatter near 
the TZ = 25.5 selection cut; and our crude astrometry, 
good to ~ 1" near the edges of survey fields. Given these 
limitations, we essentially recovered all Steidel LBGs in 
the overlapping regions. Of the 120 overlapping objects, 
all of which are cl assified photometric ally as z ~ 3 UG7Z- 
selected LBGs bv lSteidel et al.l(|2003ft . the KDF classifies 
54 as [/ ra G7?.-selected objects, 44 as BX objects, 3 as BM 
objects, and none as GTZI LBGs. Further more, of these, 
33 KD F U n GTZ LBGs have redshifts in the lSteidel et~aLl 
(2003) catalog (< z >= 2.95) as do 18 KDF BX ob- 
jects (< z >= 2. 54). The diff e rence s in classification of 
objects between ISteidel et al.1 ((2003) and the KDF can 
be attributed to photometric scatter, particularly in the 
Z7 n -band: it is reasonable that some BX objects in our 
deeper KDF catalog may have been scattered into the 
U n GTZ selection in the shallower Steidel et al. survey. 
The spectroscopic redshift distributions confirm this in- 
terpretation, as the common (KDF-selected) U n GTZ ob- 
jects are found near z ~ 3 whereas the common (KDF- 
selected) BX objects have z ~ 2.5 — a slightly higher 
reds hift than the canon ical z ~ 2.2 of the full BX sam- 
ples (ISteidel et al.ll2004ft but not surprising given the fact 
that the ([Steidel et al.1 120031 ) spectroscopy was tuned to 
target z ~ 3 galaxies rather than those at lower red- 
shifts. Furthermore, in addition to matching our objects 
with the TZum < 25.5) spectros copy of Steidel et al., we 
have obtained Gemini/GMOS dHook et al.l 12001 spec- 
troscopy for 36 UGR KDF objects, about half of them 
with TZ > 25.5. Continua are detected in all cases, al- 
though these are faint. Their faintness makes it difficult 
to readily check absorption-line spectroscopic redshifts, 
but the spectroscopy has revealed six likely LyA emitters 
below TZ = 25.5 — i.e., below the Steidel et al. survey 
limit — and another six at TZ < 25.5. No evidence of con- 
tamination by foreground interlopers was found. These 
spectroscopic tests provide additional confidence that we 
are selecting high-redshift objects by extending the Stei- 
del et al. selection criteria to fainter mag n itudes. 

As discussed in ISawicki fe Thompson! ([2006al ). at a 
rest-frame wavelength of ~1700 A, fc-corrections from 
suitably-chosen filters are minimal, and the absolute 
magnitude is very well approximated by 

M 1700 = m XobB - 5 Iog(Di/10pc) + 2.5 log(l + z). (1) 

Here, Dl is the luminosity distance and m\ obs is the 
observed magnitude in the /, TZ, GTZ, or G band for z ~ 
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4, 3, 2.2, or 1.7, corresponding to a rest-frame wavelength 
of ^1700 A. Our composite GR magnitudes are a simple 
average of the G- and 7^-band fluxes, 
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2.2. Millennium Simulation 
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In order to investigate the halo masses of galaxies sam- 
pled by our observations, we make comparisons to the 
Mille nnium Simulatior^ of dark matter structure forma- 
tion (|Springel et al.ll2005D . The simulation contains 10 10 
particles of mass 8.6 x 10 s /i _1 M Q in a cubic volume 500 
hr 1 Mpc a side, with a resolution of 5 h^ 1 kpc. The 
cosmology is a ACDM model with Hm =0.25, ft\ =0.75, 
h =0.73, and erg =0.9, similar to tha t adop ted in the 
rest of this paper. See ISpringel et all (120051 ) for a full 
description. 

We investigate the spatial distribution of dark mat- 
ter overdensities, or "halos" with Sp/p ~ 200. Data 
are available for 64 "snapshots" that were saved during 
the simulation, and from these we select all haloes with 
masses greater than ^10 n /i -1 M Q at z ~ 4, 3, 2.2 and 
1.7. 

While the majority of our results utilize the distri- 
bution of dark matter haloes in the Millennium Simu- 
lation, we also make comparisons to the predictions of 
th e semi-analytical model of galaxy formation described 
in iDe Lucia fc Blaizoti (|2007[ ) , which use the results of 
the simulation to trace the evolution of dark matter 
haloes. Model galaxy catalogs are available for the same 
64 "snapshots" as the haloes. For comparison with our 
observations, we select galaxies based on their star for- 
mation rates and total stellar masses. 

3. METHODS 

Galaxy clustering can be estimated using an angular 
two-point correlation function that measures the number 
of unique galaxy pairs at a given separation compared 
to that of a random distributio n with the same g eome - 
try. An estimator proposed by ILandv fc Szalavl (|1993l ) 
has become the standard measure of galaxy clustering at 
high-z: 

DD(8)-2DR(8) + RR(8) 

WLs{&) = rW) ' (3) 

where DD{8), DR{8), and RR(6) are the number 
of galaxy-galaxy, galaxy-random, and random-random 
pairs with angular separations between 8 — A8/2 and 
8 + A8/2. Random catalogs were compiled by inserting 
large numbers of artificial objects with various magni- 
tudes and random positions into the KDF images and 
then attempting to recover them using the same selec- 
tion procedure that is used to detect the galaxies. This 
ensures that the random catalogues have the same ge- 
ometry as the galaxy samples, avoiding areas near bright 
stars and the edges of the images. To reduce the noise 
in the random distributions, the vols calculations are 
carried out for several thousand different distributions 
drawn from the random catalogs. 



On large scales, the angular correlation function can 
be approximated by a power law of the form 



w{8) = A w 8-P. 



(4) 



Because the total number of random pairs is equal to the 
number of galaxy pairs, and wls measures the excess 
probability of finding a galaxy pair relative to finding a 
random pair, the observed correlation function cannot be 
positive for all 8. The measured correlation function is 
reduced by an amount known as the integral constraint 
(IC), which increases with clustering strength and de- 
creases with field size, 
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where the first term accounts for Poisson variance and 
the second term corrects for the fact that the mean 
galaxy density is estimated from the sample itself and no 
fluctuati ons larger than the field size have been taken into 
account ([Peebles! H980). The second term can be esti- 
mated nume rically using the random-random pair counts 
for the field (|Infantdll99l iRoche et al.lll999j ). 
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The amplitude and slope of the angular correlation 
function are then determined by the fitting function 



w{8) obs = A W 9~P - IC. (7) 

The Poissonian errors (ILandv fc Szalavlll99l are esti- 
mated for the angular correlation function, assuming the 
weak correlation limit, 



5w o bs(0) = 




(8) 



This estimate was confirmed to follow bootstrap errors 
at separations larger than 10", where the fitting is per- 
formed. 

To compare with previous studies, and to estimate the 
mass of the hosting dark matter haloes, it is useful to 
calculate the spatial correlation function, which is also a 
power law on large scales, 



(9) 



CM = - 

Jo 



where ro is the spatial correlation length and 7 = /? + 1 . 
The spatial correlation function is related to the an- 
gular correlation function through the Limber inversio n 
(jTotsuji fc Kiharalll969llMagliocchetti fc Maddoxll 1999ft . 



A 1D — 



H 7 rU F(z)rl-Uz)N 2 (z)E(z)dz 
(c/H )[fN(z)dz] 2 



(10) 



where r c (z) is the comoving radial distance, N(z) is the 
survey selection function, and 



i\ r [( 7 - 1)/2] 
r( 7 /2) 



(11) 



http: / /www. mpa-garching.mpg.de/millcnnium/ 



E{z) = yjsi m (i + z ) 3 + n A . 



(12) 
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Fig. 1. — Number-weighted selection functions for the BX, UGR, 
and GRI samples. The solid, dashed, and dash-dot lines correspond 
to magnitude ranges of TZ = 23-27, 23-25.5, and 25.5-27, respec- 
tively. 

F(z) is the evolution of clustering with redshift, which is 
assumed to be negligible within the samples considered 
here, and is set equal to 1. 

The KDF survey selection functions for different red- 
shifts and magnitude ranges are determined using the 
completeness estimates described in Sawicki & Thomp- 
son (2006). These simulations use model objects, with 
colors matching those expected among high-z galaxies, 
that are added to the data images; the same photometry 
procedure is then used on these images as on the data- 
only images to gauge the completeness as well as photo- 
metric scatter in the data. Except at z ~ 1.7 our simula- 
tions pro duce redshift d istributions that m atch very well 
with the iSteidel et ail (|1999l 120031 I2004D spectroscopic 
observations at z ~ 4, 3, and 2.2. At z ~ 1.7, the mod- 
eled selection functions do not accurately reproduce the 
observed distribution, so we have less confidence in the 
absolute values of the correlation lengths we report for 
this sample. Figure [1] shows the number- weighted selec- 
tion function estimates for various magnitude ranges for 
the four redshift samples. 

4. CLUSTERING RESULTS 

It is clear that the integral constraint has a signifi- 
cant contribution to the measured clustering strength 
of galaxies in relatively small fields. To assess the im- 
portance of the integral constraint, the correlation func- 
tion is estimated in two ways: first using the integral 
constraint calculation as outlined in the previous sec- 
tion (which we refer to as the standard IC fit), and an- 
other using an assumed value of IC=0.01 (referred to 
as the fixed IC fit), w hich is consistent with values de- 
rived from larger fields (|Adelberger et al.ll2005tlLee et al.1 
|2006t ) . Standard IC values for each sample can be found 
in Table [TJ While fixing the IC produces a somewhat 
artificial representation of the data, it can give an in- 
dication of the impact of this correction on the results 
and provide a valuable confirmation of general trends 
that may be suspect in light of possible systematics in- 
troduced by the IC. With this in mind, our analysis will 
proceed along two parallel paths, using both values of 
the IC correction. 

Since the sample sizes used in this study are generally 



too small for a robust determination of the correlation 
function slope, /?, we assume a v alue of /?=0.6, which is 
consistent with studies of LBGs (jAdelberger et ai1l2005t 
iLee et aDl2006f > and our uncertain estimates. Additional 
values of (3 = 0.8 or 0.4 are also explored for some sub- 
samples in order to to facilitate comparisons to the Mil- 
lennium Simulation. 

The left-hand panels of Figure [2] show the integral con- 
straint corrected angular correlation functions for both 
the standard and fixed IC fitting techniques for the four 
complete samples, along with the best fitting power laws. 
For the standard IC fitting procedure, the data is fit 
between 20 and 700". The lower limit minimizes the 
potential contribution of multiple galaxies sharing the 
same halo, and the upper limit roughly corresponds to 
the extent of the images. For the fixed IC fitting proce- 
dure, the lower limit is again set at 20", but the upper 
limit is reduced to ~ 100" in order to maintain a power 
law relation throughout the fitting window. Data points 
at scales larger than the outer limit are also plotted for 
completeness. Artificially reducing the IC significantly 
reduces the inferred clustering strength. The effect is 
most noticeable in samples with large angular clustering 
amplitudes (and therefore large IC values), as in the case 
of the z ~ 4 sample in Figure [5J While the IC does not 
have a large effect on the complete z ~ 2.2 and z ~ 1.7 
samples presented in Figure [2] it has a more significant 
impact on subsamples that are more strongly clustered. 

All samples show evidence of an excess on small scales 
with respect to the large scale power law, which is gener- 
ally attrib uted to multipl e galaxies residing in the same 
halo fe.g. JLee et al|[2006l ). The excess is similar to what 
has bee n observed at z ~ 4 by ILee et al.l (|2006t ) and at 
z ~ 3 bv lHildebrandt et al.1 (|2007t ). Ours is the first time 
a small scale excess has been reported for UV-selected 
galaxies at z ~ 2.2 or 1.7. 

4.1. Systematic Uncertainties 

Possibly the largest source of systematic uncertainty in 
the measured correlation functions is the IC correction. 
The relatively small field of view of the KDFs requires 
that large IC corrections need to be applied for the stan- 
dard IC fitting technique, and the resulting correlation 
lengths are sometimes significantly greater than those 
derived from larger fields. However, if a smaller IC, ap- 
propriate to a larger field is artificially implemented, the 
KDF results are in agreement with those of large area 
surveys. This means that even when the uncorrected 
measurements of two surveys are in agreement, the larger 
IC corrections for our smaller field can increase the clus- 
tering strength we infer by ~ 20 — 30%. This gives an 
indication of the potential magnitude of the systematics 
associated with the IC. It is important to note that while 
the IC can have a substantial impact on the absolute val- 
ues of the correlation lengths we derive, the important 
general trends we observe are independent of this correc- 
tion. 

Also of potential importance are the redshift selection 
functions, which are needed to perform the Limber in- 
version. In I Sawicki fc Thompson! (|2006al ). the estimated 
selection functions at z ~ 4, 3, 2.2 were shown to be 
in agreement with the Steidel et al. spectroscopic sam- 
ples, but the exact distributions, especially for the faint 
objects, remain untested. At z ~ 1.7, the simulated se- 
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Fig. 2. — a.) Angular correlation functions for the complete samples. Solid lines represent the best fitting power laws to the data between 
20 and 700" for the standard IC fitting procedure and from 20 to ~ 200" for the fixed IC fitting procedure. The data points are offset 
horizontally for clarity, b.) Angular correlation functions for various UV magnitude thresholds using the standard IC fitting technique. 



lection functions do not accurately reproduce the Stei- 
del et al. spectroscopic results. The inferred correla- 
tion lengths based on the synthetic selection functions 
are ~ 30% and ~ 20% lower for the z ~ 1.7 and z ~ 2.2 
samples and ~ 25% and ~ 20% higher for the z ~ 3 and 
z ~ 4 samples, in comparison to the correlation lengths 
based on the Steidel et al. spectroscopic results. It is also 
worth noting that the simulated selection functions show 
variations with magnitude, which cautions against us- 
ing the selection functions of the brightest, spectroscop- 
ically confirmed members, across all magnitude ranges. 
These variations generally affect the inferred correlations 
lengths by less than 10%. 

Another potential source of uncertainty is due to "cos- 
mic variance" — i.e., field-to- field variations caused by 
the fact that on small scales the universe is not homoge- 
nous. The KDF survey consists of three spatially inde- 
pendent sightlines which at least partially guards against 
the cosmic variance effect. Indeed, in our luminosity 
function work (Sawicki & Thompson 2006a) we found no 
evidence for significant field-to-field variation among the 
KDF fields. Therefore, while cosmic variance remains a 



potential concern, we do not see evidence for it in our 
data. 

Contamination of our high-redshift samples by other 
populations — most notably lower- z interlopers — could 
also affect the clustering results. In this respect it should 
be noted that spectroscopic observations at intermedi- 
ate magnitudes, i?=24-24.5 (Steidel et al. 2003, 2004), 
find only low interloper fractions of <5% in the U n GlZ- 
selected samples (z ~1.7, 2.2., 3). The interloper frac- 
tion is likely to be lower still at fainter magnitudes since 
the ratio of galaxies to Galactic stars, which are one of 
the main classes of interlopers, increases towards fainter 
magnitudes. Adelberger et al. (2005) estimate that the 
effect of interlopers is to reduce ro by up to ~ 7% in the 
BX and BM intermediate-luminosity samples and essen- 
tially by zero at z ~ 3 — a negligible effect in all three 
cases. In the G7?J-selected z ~ 4 sample the interloper 
fraction can be higher — up to 20% (Steidel et al. 1999) 
- but while this likely affects the value of ro, it is un- 
likely to affect the uncontroversial luminosity-dependent 
trend seen at this epoch. 

Finally, the choice of fitting region and bin size also 
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Fig. 3. — Correlation length, ro, with /3=0.6 for various magni- 
tude thresholds in 1Z band for the BM, BX and UGR, samples and 
in the I band for the GRI sample. The top panel is for the stan- 
dard IC fitting procedure and the bottom panel is for the fixed IC 
fitting procedure (fit between 20 and 200" for the BX and UGR 
samples and 20-150" for the GRI sample). Plus signs indicate the 
results of previous studies discussed in the text. The overall trends 
for each redshift sample are independent of the fitting technique, 
although there is a systematic offset. 



affect the clustering measurements, but these factors are 
generally found to be negligible compared to the factors 
mentioned above. 

4.2. UV Magnitude dependence 

Previous studies (e.g.. lAdelberger et al.1 120051 : 
iLee et all 120061 : iKashikawa et al.l 12006ft indicate that 
the clustering strength of UV-selected high-z galaxies 
is dependent on UV magnitude, with brighter galaxies 
clustering more strongly, but many of these studies 
have not sampled objects significantly fainter than 
M* (~ -21). The right-hand panels of Figure [H show 
the correlation functions for each redshift sample with 
different magnitude thresholds. Table [JJ and Figures |3] 
and [4] summarize the clustering parameters for the 
two IC fitting techniques. The results are presented 
for various magnitude thresholds, as well as various 
magnitude ranges at 1700 A for each sample that will 
be used to investigate the correspondence between dark 
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matter halo mass and UV magnitude in Section! 

At z ~ 3 and 4, the previous trend of increasing 
clustering strength with increasing UV magnitude (e.g.. 
Adelberger et all 120051: [Lee et al.ll200l IKashikawa et al.l 
20061) is confirmed and extended to sub-L* galaxies. At 
z ~ 2.2, this trend is only observed for objects fainter 
than 1Z ~ 25 (see TableQ] Fig. [4}. It is interesting to note 
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TABLE 1 

Correlation function parameters for UV-selected samples. 



Sample 






mag 




IC 






ro 


r a 
















[P = 0.6] 


[h~ 


L Mpc] 


[h- 1 Mpc] 


KjtXL [Z r 


A\ 

4J 


— zo.U 




T <^ Oft Q 
1 \ ZD.O 


QOQ 

ozo 


n n/in 
U.U4y 


l.UO ± U.oo 


5.2 ± 1.2 


4.4 ± U./ 


23.0 


< 


/ < 26.3 


232 


0.058 


1.25 ± 0.49 


5.5 ± 1.5 


4.7 ± 0.9 






23.0 


< 


I < 25.8 


120 


0.118 


2.59 ± 1.08 


8.1 ± 2.4 


6.4 ± 1.3 


UGR (z 


~ 3) 


23.0 




17 n 
/<- <. 27. U 


1402 


0.018 


n An _L n no 
0.40 ± U.08 


3.5 ± 0.5 


2.7 ± 0.3 


oq n 

ZO.U 




■ry ^ ofi c; 


1 nc;n 
1U0U 


n noQ 
U.Uzo 


n ftQ n 11 
U.OO ± U.ll 


4.5 ± 0.5 


no 1 no 
O.Z ± U.O 






23.0 


< 


7? < 26.0 


709 


0.039 


0.73 ± 0.16 


5.4 ± 0.6 


3.7 ± 0.4 






23.0 


< 


TL < 25.5 


424 


0.055 


1.23 ± 0.27 


6.4 ± 0.9 


4.3 ± 0.7 


BX (z ~ 


2 2) 


23.0 


< 


TL < 27.0 


2312 


0.015 


0.35 ± 0.05 


3.6 


± 0.3 


3.1 ± 0.1 


23.0 


< 


ft < 26.5 


1910 


0.028 


0.65 ± 0.06 


5.1 ± 0.3 


3.9 ± 0.2 






23.0 


< 


TL < 26.0 


1483 


0.037 


0.84 ± 0.08 


5.8 


± 0.4 


4.1 ± 0.2 






23.0 


< 


•R. < 25.5 


1020 


0.038 


0.86 ± O.ll 


5.8 


± 0.5 


4.1 ± 0.3 


BM (z - 


' 1-7) 


23.0 


< 


TL < 27.0 


1966 


0.026 


0.59 ± 0.05 


3.6 


± 0.2 


3.0 ± 0.1 


23.0 


< 


TL < 26.5 


1707 


0.033 


0.76 ± 0.06 


3.9 ± 0.2 


3.1 ± 0.1 






23.0 


< 


7£ < 26.0 


1330 


0.034 


0.76 ± 0.08 


3.5 ± 0.2 


2.8 ± 0.2 






23.0 


< 


TL < 25.5 


971 


0.035 


0.79 ± O.ll 


3.2 ± 0.3 


2.6 ± 0.2 


GRI (z r 


-4) 


23.0 




/ < 25.8 


120 


0.118 


Z.49 ± l.Uo 


8.1 


± 2.4 


6.4 ± 1.3 


OK Q 
zo.o 




r . ofi q 

i \ ZD.O 


1 TO 

i i y 


n nQO 
U.Uoz 


I. ( ( ± U.D4 


7.1 


± 1.7 


K A J- 1 1 
0.4 ± 1.1 






25.8 


< 


/ < 26.8 


201 


0.056 


1.18 ± 0.56 


6.0 


± 2.0 


4.8 ± 1.2 


UGR (z 


~ 3) 


ZO.U 




/V 2_ ZO.O 


/(OK 

4Z0 


U.UOO 


i oq _l_ n no 
l.zo ± U.uy 


6.4 


± 0.9 


An 1 r. rr 
4.0 ± U. 1 


ok n 
ZO.U 




"7T> Oft n 


K1 Q 

Olo 


fl OQ7 

U.Uo / 


n qo _l_ n oo 
U.oz ± U.zz 


5.3 


± 0.9 


q p. A- n p. 
o.D ± U.D 






or; n 
ZO.U 


< 


n> << 07 n 
/<- v. z/.U 


1ZUU 


U.UZ4 


n kq _L n o*7 
U.Oo ± U.z / 


1.3 


± 0.5 


on 1 no 

z.y ± U.O 


BX (z ~ 


2.2) 


23.0 


<. 


G /<. <- zo.U 


561 


0.037 


n o q _L n o 1 
0.8o ± U.21 


5.5 


± 0.9 


4.7 ± 0.6 


23.0 




Lr !<, <- ZO.O 


943 


0.042 


n o k j_ n 1 o 
0.90 ± U.lz 


6.1 


± 0.5 


4.7 i 0.3 






24.0 


< 




838 


0.052 


1.19 ± 0.14 


7.1 


± 0.5 


5.3 ± 0.4 






25.0 


< 


Gfc < 26.0 


833 


0.052 


1.19 ± 0.14 


7.6 


± 0.6 


5.6 ± 0.4 






25.5 


< 


GR < 26.5 


906 


0.025 


0.55 ± 0.13 


5.0 


± 0.7 


4.5 ± 0.4 






26.0 


< 


GR < 27.0 


874 


0.012 


0.25 ± 0.13 


3.3 


± 1.2 


3.3 ± 0.6 


BM (z - 


'1-7) 


23.0 


< 


G < 25.0 


486 


0.020 


0.42 ± 0.25 


2.0 


± 0.9 


1.7 ± 0.5 


23.5 


< 


G < 25.5 


795 


0.017 


0.37 ± 0.14 


2.0 


± 0.5 


1.9 ± 0.3 






24.5 


< 


G < 25.5 


576 


0.013 


0.26 ± 0.16 


1.8 


± 0.8 


1.8 ± 0.6 






25.0 


< 


G < 26.0 


703 


0.019 


0.41 ± 0.16 


2.8 


± 0.7 


2.5 ± 0.5 






25.0 


< 


G < 27.0 


1390 


0.025 


0.56 ± 0.08 


1.1 


± 0.4 


3.6 ± 0.3 






25.5 


< 


G < 26.5 


744 


0.032 


0.72 ± 0.16 


5.0 


± 0.7 


4.6 ± 0.4 






26.0 


< 


G < 27.0 


683 


0.032 


0.72 ± 0.16 


6.0 


± 0.9 


5.1 ± 0.5 



a Correlation length for fixed IC value of 0.01. 

that the most clustered objects at z ~ 2.2 are not the 
brightest ones, but those with K = 25 -26 (Fig. H. This 
is in d isagreement with the findings of lAdelberger et aLI 
(2005) for objects with 1Z < 25 but i t is in agreement 
with the findings of Q uadri et al.l ()2007[ ). who report that 
for their sample of if-selected galaxies at z ~ 2.5, those 
with R > 25 cluster more strongly than those with R < 
25. Our inclusion of fainter optically selected objects re- 
veals an unexpected similarity in the clustering proper- 
ties of if -selected galaxies and optically selected galaxies 
at z ~ 2.2. At z ~ 1.7, the UV magnitude trend ap- 
pears to be reversed with respect to z ~ 3 and 4, with 
the faintest galaxies being the most strongly clustered. 
It should be noted that while the numerical values of ro 
vary, the overall trends in each sample are independent 
of the IC fitting technique. We discuss the potential im- 
plications of these trends with redshift in the context of 
galaxy downsizing in Section 16.11 

4.3. UV color dependence 

One relationship that has not been previously reported 
is an association between the correlation length of a sam- 
ple and its UV color. Figure [5] shows the relation be- 



tween sample UV color threshold and correlation length 
for galaxies with M1700 < —19. At z ~ 4, 3, and 2.2, 
the galaxies with blue UV colors are significantly more 
clustered. There does not seem to be a relation among 
UV-bright galaxies at z ~ 1.7, but faint galaxies with UV 
red colors show evidence of being more strongly clus- 
tered. This apparent UV color-dependent clustering is 
independent of the IC fitting technique, and provides a 
caution against making direct comparisons between stud- 
ies that use different selection techniques, or between 
galaxy samples at different rcdshifts, where the exact 
color selections are not identical. We note that finding 
redder galaxies to be weakly clustered seems counterin- 
tuitive, and we discuss the possibility that thi s tre nd is 
due to differences in dust properties in Section [6.21 

The especially striking UV color relationship at z ~ 
2.2, in which the measured correlation lengths for both 
IC fitting techniques increase to well in excess of 10 h~ l 
Mpc (for the bluest samples) may seem unusual for UV- 
selected galaxies, but we are reminded of the similarities 
mentioned in the previous section to if -selected galaxies, 
which are routinely observed to have correlation lengths 
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of these magnitudes (jQuadri et al.l2007l) . Taken together 
with the observed magnitude dependence of the cluster- 
ing, this provides evidence of significant overlap between 
z ~ 2.2 optically selected samples and those obtained 
using if-selection, which are generally thought to be 
very different populations due to differences in clustering 
strength. 

4.4. Comparison with Other Surveys 

There are still discrepancies between various surveys 
concerning the clustering properties of high redshift 
galaxies, which are generally attributed to differences in 
selection techniques. This suggestion is supported by the 
evidence presented in Section 14.31 that indicates that in 
addition to the UV magnitude of objects, UV color can 
also significantly affect the observed clustering proper- 
ties. 

One advantage of the KDF survey is that it shares the 
same selection criteria as work by the Stcidcl et al. group, 
a fact that allows us to make comparisons that are likely 
independent of selection differe nces. For the same selec - 
tion criteria as those used here, lAdelberger et al.l (|2005l ) 
report correlation lengths of 4.2 ± 0.5 and 4.0 ± 0.6 
Mpc for galaxies at z ~ 2.2 and 3 with 7 ~ f.6 and 
1Z = 23.5-25.5. Despite the fact that we have identical 
selection techniques, when the KDF data of similar mag- 
nitudes are fit with the standard IC fitting procedure, we 
obtain larger results 5.8 ± 0.5 and 6.4 ± 0.9 Mpc. 
The discrepancy can be reduced if force smaller IC val- 
ues in our analysis. If the fixed IC fitting technique is 
used, the lower IC value reduces the correlation lengths 
to 4.1 ± 0.3 a nd 4.3 ± 0.7 h -1 Mpc, w hich are consistent 
with those of lAdelberger et all (|2005l) . This emphasizes 
the importance of the IC. Also worth mentioning is that 
for galaxies with t he sam e m agnitude range at z ~ 3, 
iHildebrandt et al.l (|2007| ) and lYoshida et all ([20081 ) re- 
port values of 5.0 ± 0. 2 and 5.5 ± 0.9 fe" 1 Mp c , which 
are larger than those of lAdelberger et all (2005), but still 
slightly less than our stan dard IC fit t ing re sults. 

For the GOODS fields, ILee et al.l ([200(1 report clus- 
tering measurements for galaxies at z ~ 3 and 4. At z ~ 
3 they report a value of rp=4.0 ± 0.2 h~ x Mpc, similar 
to the value reported in lAdelberger et al.l ([2005) . Our 
larger results for the standard IC fitting proc edure can 
again be traced to large IC values. At z ~ 4, ILee et all 
(2006) obtained a deep sample with z 8 5o < 27, and re- 
port a correlation length of 2.9 ± 0.7 h^ 1 Mpc. For a 
similar magnitude cut (/ < 26.8), we obtain a value of 
5.2 ± 1.2 h^ 1 Mpc with the standard IC fitting proce- 
dure and a value of 4.4 ± 1.2 h^ 1 Mpc for the fixed 
IC fitting procedure. This difference could possibly be 
caused by a difference in sample selection. In addition 
to being magnitude selected in different filters, it is also 
possible that the Lee et al. sample contains redder ob- 
jects which are shown to cluster much less strongly in 
Figure [5j In the KDFs, objects in the z ~ 4 sample 
with 1Z — I > 0.1 are essentially not clustered, and they 
significantly reduce the correlation length for the total 
sample. Another difference between the samples is that 
at the limit of Zgso < 27, the GOODS data are estimated 
to be ~10% complete, while the KDFs are estimated to 
be ~50% complete at I < 26.8. 

The KDF results at z ~ 3 and 4 are in good agree- 
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IC fitting procedures. Histograms indicate the cumulative number 
of galaxies below the color threshold and M1700 < —19- The open 
points in the z ~ 1.7 panel are for all galaxies with —20 < Mi 700 < 
— 19 redder than a given threshold in G — TZ. 



ment with thos e of t he Subaru Deep Field (SDF). 
Kashikaw a et al.l (|2006D obtained az~4 sample with a 
limiting magnitude of «=27.43. They report single power 
law fits to measurements on all scales, including r < 20". 
If our fitting method is altered to include all scales, and 
we assume the same values of /3 that they report, the 
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KDF results for the standard fitting technique agree well 
with their findings. As previously mentioned, the z ~ 3 
resul ts from the SDF, wit h a limiting magnitude of z' = 
25.5 (|Yoshida et alJl2008l ). also agree with our results. 

The correlation lengths repor ted bvlOuchi et al.1 (12005 ) 
for z ~ 4 objects, and those of iHildebrandt et al.l (2009) 
at z ~ 3 and 4 are generally consistent with our results 
for the fixed IC fitting procedure, but are smaller than 
our results for the standard IC fitting procedure. 

In summary, the KDF clustering measurements are 
generally found to be in good agreement with previous 
studies, allowing for differences in selction and fitting 
techniques. 

5. COMPARISON WITH MILLENNIUM SIMULATION 

The spatial correlation function of dark matter halos is 
known to be a function of halo mass (e.g., IMo & Whiti 
1996), and it has also been established that the corre- 
lation functions of UV -selected galaxies are dependent 
on UV magnitude (e.g., Adclbcrg er et al.f 2005; Lee et al.l 
2006ir In order to investigate a 



100 



2006; Kashik awa et al 



possible relationship between UV magnitude and halo 
mass, we use the publicly- available Millennium Simu- 
lation halo catalogs to estimate the spatial correlation 
function of halos of various mass, and relate them to our 
observations of UV-selected galaxies. 

The dependence of halo clustering on mass in the Mil- 
lennium simulation is investigated in two manners. In 
one method, all haloes with masses greater than a given 
threshold are selected, and for the second method, all 
haloes within a given mass range are selected. The first 
method enables a comparison between the number den- 
sities of haloes and observed galaxies, and the second is 
used to develop a relationship between UV magnitude 
and a typical halo mass. The results of the two tech- 
niques are illustrated in Figure |6] and Figure [7j respec- 
tively, and compared to the observed clustering of galax- 
ies with various UV magnitude thresholds and ranges. 
To exclude the effects of multiple galaxies residing in 
the same halo, the comparison is made on large scales 
(1 < r < 10 h^ 1 Mpc). Given the large uncertain- 
ties, rough halo mass estimates are established by eye 
for both techniques, rather than implementing a formal 
fitting procedure. In cases where the assumed slope of 
/3=0.6 for the observed galaxies does not agree with the 
predictions of the simulation, a slope of 0.4 or 0.8 is also 
investigated. The quoted uncertainties for all mass esti- 
mates are determined by plotting the associated correla- 
tion function uncertainties over the simulation results in 
the same manner as in Figures [6] and [7j 

5.1. Halo Occupation Numbers 

One advantage of comparing our clustering measure- 
ments to the distribution of dark matter haloes in the 
Millennium Simulation is that it allows us to make com- 
parisons between the observed number densities of galax- 
ies and the number densities of haloes with similar cor- 
relation functions. The comparison is characterized by 
the mean halo occupation number, N occ , which is de- 
fined as the ratio of the galaxy number density to that of 
the haloes. The number densities of galaxies sampled by 
the KDFs are found by inte grating the observed luminos- 
ity functions presented in iSawicki fc Thompson! (|2006al) 




r (h 1 Mpc) 

Fig. 6. — Comparison of observed spatial correlation lengths 
for various UV magnitude thresholds (summarized in Table [TJ 
to the predictions of the Millennium Simulation for different halo 
mass thresholds. Filled points correspond to increasing halo mass 
thresholds of Iog(fc- 1 M )=1O.7, 10.9, 11.1, 11.3, 11.5, 11.7, 11.9, 
12.1. Solid lines represent the results of the standard IC fitting pro- 
cedure, and the dashed lines are for the fixed IC fitting procedure. 



down to a given UV magnitude threshold. The corre- 
lation functions for these galaxies are then compared to 
the predictions of the Millennium simulation for various 
threshold masses (Figure [5]), which have known number 
densities. Table [5] summarizes the number densities of 
observed galaxies and dark matter halos with similar cor- 
relation functions. The ranges reported for the quantity 
nobs/nhaio are found from the maximum and minimum 
values allowed by the uncertainty in n^ Q ; . 

For most samples, the number of estimated haloes is 
generally consistent with the number of galaxies, imply- 
ing that there is about one galaxy for every halo. This is 
not the case for the z ~ 2.2 and z ~ 1.7 samples (espe- 



10 



Savoy et al. 



100 



10 



0.1 



E 

- • 


1 1 


1 


1 


i 111 = 

GRI : 






















• - 
• 




• 






• 






• 






• • 
• 












• 


> 


----- »-V"- : C"":- 








• 


• 










• 



10 



1 r 



0.1 



10 



0.1 



10 



1 r 



0.1 




r (h 1 Mpc) 

Fig. 7. — Comparison of spatial correlation lengths of UV lumi- 
nosity ranges (summarized in Table [TJ to the predictions of the 
Millennium Simulation for distinct halo mass ranges. Filled points 
correspond to increasing halo mass ranges of log(/i — 1 Mq)=10.7- 
11.0, 11.0-11.3, 11.3-11.5, 11.5-11.7, 11.7-11.9, 11.9-12.1, 12.1-12.3, 
12.3-12.5. An additional mass range of log(h — 1 Mq)=12.5-12.7 is 
shown for the BX and BM samples. Solid lines represent the re- 
sults of the standard IC fitting procedure, and the dashed lines are 
for the fixed IC fitting procedure. 



daily for the standard IC fitting results) , where there are 
more galaxies than haloes by a factor of a few. We are 
again reminded o f possible simi l aritie s to th e if-selected 
z ~ 2.2 sample of lQuadri etafl (|2007l l200l . who report 
a very large halo occupation number of ^40 for their red- 
dest sample (DRGs, with J — K > 2.3), and a value of 
~l-4 for their complete sample (K < 21). 

5.2. Relation Between Halo Mass and SFR 

Comparing the observed correlation functions of ob- 
jects with distinct UV magnitude ranges to the predicted 
clustering properties of distinct halo masses (Figure [7]) 
allows us to investigate UV luminosities and SFRs asso- 



ciated with halos of a given mass, as a function of red- 
shift. The left-hand panels of Figure [8] show, for the 
first time, the estimated relationships for galaxies at z ~ 
1.7, 2.2, 3, and 4 for the standard and fixed IC fitting 
procedures (neglecting dust absorption). It can be seen 
that for both IC fitting techniques, at z ^3 and 4 the 
dark matter mass-UV magnitude relationships are es- 
sentially indistinguishable, but at z ~ 2.2 the relation 
is significantly below the higher redshift results. This 
could possibly indicate that star formation was more ac- 
tive at higher redshift in halos o f a given mass, which 
is consistent with the findings of iLee et al.l ((2006) , who 
find that star formation was more efficient at z ~ 5 than 
at z ~ 3-4. It is worth noting that these conclusions are 
subject to the amount of dust absorption in the galax- 
ies, but they are independent of the IC fitting procedure 
used. 

The predicted relationship between halo mass and 
UV lumin osity from the Millennium Simulation galaxy 
catalogs (|De Lucia fc BlaizoD [2007) was also investi- 
gated. To make the comparison, the predicted star 
formation rates of the modeled galaxies were first con- 
verted to M1700 magnitudes u sing the following relation 
(jSawicki fc Thompson|[2006al) : 



SFR = 6.1 x l0-( 8+0 - 4Miro °)M ( 7 



(13) 



which neglects dust absor ption and is ap propriate for 
a stellar population with a Sal peterl (|1955l ) initial mass 
function over a range of 0.1 < M/Mq < 100. 

In order to make comparisons between the simulated 
galaxies and our results, at each redshift the predictions 
are reduced to a median relation between halo mass and 
UV magnitude (represented as solid lines in Figure HJ . 
While the median values are convenient for making com- 
parisons, it should be noted that the range of halo masses 
are generally symmetrical in log(MoM) and typically 
have a standard deviation of a factor of ^2 in mass. 

The absence of dust in the conversion is clearly inac- 
curate, and as expected, the predicted M1700 magnitude 
for galaxies of a given halo mass is much greater than 
our observed relations. However, magnitude-dependent 
dust corrections are now available at z ~ 2.2 from U- 
through-i/ SED fitt ing of UV-selected galaxies in th e 
Hubble Deep Field (jSawicki et al.l [20071: Ifca^Ickll I20TI . 
and we apply these to our presents KDF results. It is as- 
sumed that the magnitude-dependent corrections at the 
other redshifts are identical to those at z ~ 2.2. 

Accounting for the absorption by dust significantly al- 
ters the observed relations (right-hand panels of Figured]) 
and brings our results into much better agreement with 
the predictions. After the corrections are made, the z ~ 
3 and 4 results for both IC fitting procedures agree within 
the errors with the predictions of the Millennium Simu- 
lation. The agreement is especially striking, given the 
considerable uncertainties in both the models and our 
work. At lower redshifts there are large differences be- 
tween our results and the model predictions, perhaps 
reflecting changes in the evolution of galaxies at these 
redshifts that have not been considered in the model. 

5.3. Stellar to Dark Matter Mass Relation at z ~ 2.2 

Another interesting quantity that can be investigated 
is the relation between stellar and dark matter halo mass, 
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TABLE 2 
Galaxy and halo number densities. 







^obs 


Mass 




Mass 


nhalo 


Sample 


mag 


(xicr 3 ) 


Threshold a 


(xl(T 3 ) n obs /n halo a 


Threshold 15 


(xl0~ 3 ) n obs /n ha i h 






[Mpc" 3 ] 


[logCh-iMe)] 


[Mpc" 3 ] 


[logCh-iMe)] 


[Mpc" 3 ] 



GRI 23 < / < 26.8 1.5 

23 < / < 26.3 1.0 

23 < / < 25.8 0.6 

UGR 23 < TZ < 27.0 4.8 

23 < TZ < 26.5 3.2 

23 < TZ < 26.0 2.0 

23 < TZ < 25.5 1.2 

BX 23 • 7C 27.0 7.5 

23 < TZ < 26.5 5.6 

23 < TZ < 26.0 3.9 

23 < TZ < 25.5 2.6 

BM 23 • 7C 27.0 24.1 
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Fig. 8. — Relation between dark matter halo mass and UV luminosity/SFR (left and right axes) for all redshift samples. The bottom 
panels correspond to the standard IC fitting procedure and the top panels to the fixed IC fitting procedure. The panels on the left are 
the observed relations, and the panels on the right are the results after corre ctions for dust attenuation have been made. The solid lines 
correspond to the predictions from the Millennium Simulation galaxy models (Dc Lucia & Blaizot 2007), reduced to median relations. 
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which can be estimated by combining the above relation 
between halo mass and UV luminosity with a recently 
developed relation between stellar mass and UV lumi- 
nosity f or star forming galax ies at z ~ 2.2 ()Sawickill2011l : 
see also Sawicki ct al. 2007): 
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log- 
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700, 
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where Mi 700 is the absolute 1700 A luminosity not cor- 
rected for dust. This relation is derived from broad- 
band spectral energy distribution (SED) fitting for a 
sample of UV-selected z ~ 2.2 galaxies extending 
to i?,=28 in the Hubble Deep Field. The observed 
^300^450^606^814^110^160 S EDs are compared t o the 
spectral synthesis models of iBruzual fe Charlotl (|2003f ) 
with a constant star formation rate and a metallicity 
of 0.2 Zp,. The s ynthe tic spectra are reddened using 
the lCalzetti et all ( 20001 ) starburst extinction curve, and 
shifted in wavelength and intensity according to the 
adopted cosmological parameters, which are the same 
as those used i n this study. The relation is valid for the 
iSalpeterl (jl955[ ) stellar initial mass function (IMF). 

Combining the z ~ 2.2 points in Figure [8] with Equa- 
tion [14] produces the results shown in Figure |U The two 
sets of points correspond to the two IC fitting procedures, 
and the greyscale represents the pr edictions of the Mil- 
lenniu m Simulation galaxy models of lDe Lucia fc Blaizotl 
(|2007h . There is a significant disagreement, particularly 
at intermediate halo masses, between the results for both 
fitting procedures and the models and our stellar mass 
estimates are low by a factor of about 5-10, compared 
to the models. This discrepancy could be indicative of 
limitations in: (1) our method of estimating dark mat- 
ter mass from clustering measurements, (2) the SED fit- 
ting process used to estimate the stellar masses, or (3) 
the galaxy models. One possible explanation for the dis- 
crepancy is that the total stellar mass is underestimated 
because it is derived exclusively from observations of 
the UV continuum which only sample young, recently- 
formed stars, and do not account for the contribution 
of existing older stellar populations. Another possibility 
is that the estimated dark matter masses are biased to- 
wards massive haloes that strongly affect the clustering 
measurements, while the stellar mass estimates are bi- 
ased towards the more numerous and less-massive, but 
still actively star -forming galaxies. We note, however, 
that applying Eg . 1141 to the clustering measurements of 
lAdelberger et al.l (|2005D gives a reasonably good agree- 
ment with the models, albeit this can be tested only at 
the highest halo masses. 

6. DISCUSSION 

6.1. Magnitude-Dependent Clustering and Galaxy 
Downsizing 

The premise of galaxy downsizing originated with 
the observation that the sites of highest star formation 
shift from high to low stellar mass systems with time 
(jCowie et al.| [l996). A logical extension of this scenario 
would seem to be that star formation first began in the 
most massive dark matter haloes and then, as time pro- 
gressed, the sites of highest star formation migrated from 
higher to lower mass haloes. Such an extension of galaxy 
downsizing from stellar to halo mass is suggested by 
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Fig. 9. — Relation between stellar mass and dark matter halo 
mass at z ~ 2.2 for the different IC fitting techniques. The 
greyscale represents the predictions from the Millennium Simula- 
tion gala xy catalogs. For compari son, the green triangles show the 
results of Adelbcrgcr ct al. (2005). 



IHeinis et al.l (|2007|) , who use clustering measurements of 
UV-bright star forming galaxies at z ~ 1 and 0.2 to deter- 
mine that these objects reside in significantly lower mass 
haloes than at z > 2. Our magnitude-dependent cluster- 
ing results at z ~ 4, 3, 2.2, and 1.7, are well-suited to 
extend the framework of galaxy downsizing to halo mass. 

At z ~ 3 and 4, the UV-bright galaxies in our sam- 
ple are observed to cluster more strongly than UV-faint 
ones, in agreement with previous studies. We propose 
that at these redshifts, the current UV-limited samples 
consist mainly of massive dark matter haloes, in which 
there is a strong relation between halo mass and SFR 
(as indicated by UV magnitude). By z ~ 2.2, a differ- 
ent picture begins to emerge, as we find evidence that 
the most strongly clustered galaxies in our sample are 
no longer those with the highest UV luminosities. By 
z ~ 1.7, the UV magnitude relation appears to be com- 
pletely reversed, with the UV-faint galaxies in our sample 
clustering more strongly than UV-bright ones (Fig. @|. 

Our results suggest that for the objects in our sample, 
the most massive dark matter haloes host galaxies with 
high SFRs (M1700 < -20, or > 50 M yr" 1 ) at z ~ 3 
and 4, moderate SFRs (-20 < A/1700 < -19, or - 20 M Q 
y^ 1 ) at z ~ 2.2, and lower SFRs (-19 < M 1700 < -18, 
or ~ 2 M© yr -1 ) at z ~ 1.7. The decrease in star for- 
mation we observe in the most m assive haloes a t z ~ 2 
is c onsistent wit h the findings of lHartlev et al.l (|2008l ) 
and iBlanc et~al1 ([2008) who report that the popula- 
tion of passively evolving BzK-selected galaxies is more 
strongly clustered than are star-forming galaxies at this 
redshift. In this scenario, by z ~ 2, the epoch of large- 
scale star formation is coming to an end in many of the 
most massive haloes, which is reflected in the larger cor- 
relation lengths of UV-faint galaxies in our sample. Our 
results are also consistent with the weak relation be- 
tw een UV magnitude and clustering strength observed 
by IHeinis et al.l (|2007D at z ~ 1, because our results sug- 
gest that the absolute UV magnitude of the galaxies in 
the most massive haloes should be Mnoo > —18, below 
the threshold of their survey. The implied shutdown of 
star formation in the most massiv e haloes at 1 < z < 2 i s 
also consistent with the findings of lArnouts et al.1 (2007), 
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who report a buildup of the quiescent red sequence of 
galaxies between 1 < z < 2. The developing consen- 
sus that z ~ 2 is the redshift at which star formation 
is being quenched in the most massive haloes, while 
lower mass haloes continue to rapidly form new stars, 
is further supported by our magnitude-dependent clus- 
tering measurements. These results are also supported 
by the UV luminosity and SFR dens ities derived for the 
KDFs (|Sawicki fc Thompson! I2006bf ). which show that 
the dominant contributor to the total luminosity changes 
with redshift, with less luminous galaxies becoming more 
dominant as time progresses. 

It has occurred to us that the number densities of the 
strongly clustered faint objects in the z ~ 1.7 sample 
are much larger than those of other strongly clustered 
populations, which is a concern if they are all believed 
to represent the most massive haloes. However, we have 
used results from the Millennium Simulation to deter- 
mine that the observed ~ 50% increase in clustering for 
faint objects can be attainted by adding about 1 in 5 
high-mass haloes to a sample of weakly clustered low- 
mass haloes. It seems plausible that high-mass haloes 
could account for ~ 20% of the faint objects in our sam- 
ple, which could raise the clustering strength by the ob- 
served ~ 50%. It should be noted that these estimated 
increases in clustering strength only apply to very weakly 
clustered "background" populations and should not be 
generalized to other populations. 

We believe our clustering measurements provide a new 
line of evidence for galaxy downsizing. By using cluster- 
ing measurements to determine the typical dark matter 
halo mass of star forming galaxies across a broad range 
of luminosities and redshifts, we can trace the evolution 
of the sites of highest star formation (as indicated by 
UV luminosity) from high-mass haloes at high redshift, 
to lower-mass haloes at lower redshift. This suggests an 
observational link between galaxy downsizing and halo 
mass. 

The reversal of the UV magnitude trend at z ~ 2 
has not been previously noted due to the lack of deep 
UV coverage needed to select galaxies at these interme- 
diate redshifts. The very deep coverage of the KDFs 
across wavelengths needed to efficiently select galaxies 
from z ~ 1.5 — 4.5, has allowed us to detect this impor- 
tant trend in the evolution of galaxies in the context of 
halo mass. The beginnings of this trend were first hinted 
at by iQuadri et al.l (j2007l ). who report that objects in 
their iC-selected sample at z ~ 2.5 with R > 25 cluster 
more strongly tha n those with R < 2 5. The concept is 
also supported bv iHeinis et al.l (|2007l) . who report that 
there is very little relation between UV magnitude and 
clustering strength for galaxies at lower z. Our findings 
provide a significant confirmation and extension of these 
results. 

6.2. UV color relation: Observing the 2175 A UV 
absorption bump? 

In Section l4~3"l we report that galaxies with blue UV col- 
ors (in G— TZ or TZ— I) appear to cluster more strongly at 
z ~ 4, 3, 2.2, with the relationship being most striking at 
z ~ 2.2. This relationship does not exist at z ~ 1.7 and 
could even be reversed, with UV-red galaxies possibly 
clustering more strongly (Figure [5]). We have verified 
that this trend also exists in BzK -selected galaxies at 



z ~ 2 by using the publicly available MUSYC catalogs 
for the ECDFS field, to confirm that galaxies in that 
sample with blue UV colors (in B — R) also cluster much 
more strongly (but not for R — I). The trend of objects 
with blue UV colors clustering more strongly is opposite 
to what is observed in optical colors. IQuadri et al.1 ()2007f ) 
report that for their i^-limited sample at z ~ 2.5, opti- 
cally red galaxies (in J — K) cluster much more strongly, 
and raise the question of whether this effect is due to old 
stellar populations or dust. 

An intriguing possibility is that the strong relation we 
observe between UV color and clustering strength is a re- 
sult of the well-known 2175 A UV absorption bump. Lo- 
cally, this broad spectral feature is prominent along most 
sight-lines in the Milky Way galaxy, is moderately strong 
in the Large Magellanic Cloud (LM C), and is absent in 
the Small Magellanic Cloud (SMC) (IGordon et al.ll200l 
(see Figure [T0|) . The f eature has not been detected in 
local starburst galaxies ( Calzctt i et al.lll~994D . and in the 
Milky Way it has been associated with regions of diffuse 
ISM a nd seems to be sup pressed in dense star-forming 
areas (Whiffet et al. 2004). One of the currently favored 
carriers for the absorption band is polycyclic aromatic 
hydrocarbons (PAHs) , which are produced by AGB stars 
of intermediate ages. An age dependence for this absorp- 
tion feature could explain its absence in areas of recent 
star formation. 

At z ~ 2.2, this distinct absorption feature has recently 
been observed in composite spectra of UV-bright galax- 
ies with redder overall UV continuum slopes, b ut not 
in those with blue UV slop es dNoll fc Pierinil l 2005l ). This 
work has been extended bv lNoll et al.1 (|2007l ). who report 
that at z ~ 2.2, the feature is stronger in more "mature" 
galaxies, consistent with the carrier of the absorption be- 
ing produced by older stellar populations. 

The argument for this feature playing a role in our 
UV color relation is especially convincing for the BX 
sample at z ~ 2.2, where the rest-frame wavelength of 
2175 A falls directly in the TZ band (see Figure [ID]). If 
a strong absorption bump is present, it will significantly 
reduce the TZ band flux and actually make the observed 
G — TZ color bluer compared to galaxies without the fea- 
ture. Since the absorption bump seems to be associated 
with relatively evolved stellar populations, it follows that 
it would be more prominent in more massive haloes at 
z ~ 2.2, in which star formation is shutting down. In 
this scenario, galaxies with relatively blue G — TZ colors 
would then be the most massive and strongly clustered, 
in agreement with our observations. 

In studies of high-z galaxies, it is usually as- 
sumed that the dust properties are si milar to those 
found in local starburst galax ies (e.g.. iSawicki fc Yeel 
(19981 : iPapovich et alJ (|200l : iShaplev et al.l (|2001l 
120051 ) : lYabe et al.l (|2009| )). Since it seems possible that 
old high-z galaxies contain a Milky Way-type mixture of 
dust, it will be interesting to consider what effects this 
could have on global properties, such as SFRs, that are 
derived from SED-fitting. 

The redshift dependence of the UV-color relations we 
observe can actually support the concept of galaxy down- 
sizing discussed in the previous section if we believe that 
the 2175 A absorption bump is driving the trends and 
that it is more prominent in massive haloes. If UV color 
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Fig. 10. — Average extinction curves for local galaxies compared 
to our UnGTZI filter set at different sample redshifts. Note the 
prominence of the 2175 A absorption bump in Milky Way dust, 
and its locatio n in the 1Z band at z ~ 2. The data for the LMC and 
SMC are fr om IGordon et al.l 12003) , the MW curve (with R v = 
3.1) is from Cardclli et al. (1989), and the starburst curve (with 
R v = 4.05) is from Calzctt ietal] J2000f) . 



is an indicator of halo mass (which seems possible, given 
the strong clustering dependence shown in Figure [5]), our 
observations can be interpreted in the following way: (1) 
at z ~ 3 and 4, there are few low mass haloes in the 
sample, so the UV color dependence is relatively small, 
but can still be used to select the most massive haloes. 
(2) At z ~ 2.2, the sites of highest star formation are 
beginning to migrate from high to low mass haloes, re- 
sulting in a significant range of halo masses in the sample, 
and the ability to separate them based on color produces 
a very large UV-color relation. (3) At z ~ 1.7, most 
galaxies in the sample are hosted by low mass haloes in 



which the absorption bump is not prominent, and the 
UV-color relation is removed or reversed, with UV-red 
galaxies possibly clustering more strongly. This scenario 
provides additional support for galaxy downsizing. 

Our proposal that the UV-bright objects at z ~ 2 are 
hosted in a mixture of high and low mas s haloes is consis- 
tent with results of iConrov et al.l (|2008T ). who argue that 
based on the most up-to-date number density observa- 
tions, star forming galaxies at z ~ 2 (specifically, those 
with 1Z < 25) cannot evolve into massive red galaxies, a s 
has been previously suggested (jAdelberger et all 12005). 
Taking into account the merging history of z ~ 2 star 
forming galaxies, they find that the number densities are 
most consistent with those of typical ~ L* galaxies at 
z ~ 1 and 0. Our interpretation allows for a fraction of 
the z ~ 2 star forming galaxies (those with a prominent 
2175 A absorption bump and blue G — R colors) to re- 
side in the most massive haloes and possibly evolve into 
ellipticals by z ~ 0, while the majority may evolve into 
typical L* galaxies like our own. 

7. SUMMARY AND CONCLUSIONS 

We have investigated the clustering properties of UV- 
selected galaxies at z ~ 4, 3, 2.2, and 1.7, extending to 
those significantly below L*. The main findings of our 
study can be summarized as follows: 

1. At z ~ 4 and 3, we find that the UV-brightest 
galaxies are the most strongly clustered, in agree- 
ment with previous studies, but this is not the case 
at z ~ 2.2 and 1.7. We believe that we are ob- 
serving galaxy downsizing, as the sites of highest 
star formation migrate from high mass dark mat- 
ter haloes at high redshift to lower mass haloes at 
lower redshift. 

2. By comparing our clustering measurements to the 
predictions of the Millennium Simulation, we de- 
rive halo mass estimates for the galaxies in our 
samples, and from these comparisons, we develop 
relationships between halo mass and SFR for each 
sample. Our results are consistent with a shut- 
down of star formation in massive haloes as time 
progresses. In particular, we estimate that the 
the most massive dark matter haloes in our sam- 
ple (~ 1 — 5 x 10 12 M ) host galaxies with high 
SFRs (Mi 700 < -20, or > 50 M Q yr" 1 ) at z ~ 
3 and 4, moderate SFRs (-20 < A/1700 < -19, 
or ~ 20 M yr" 1 ) at z ~ 2.2, and lower SFRs 
(-19 < M 1700 < -18, or - 2 M yr" 1 ) at z ~ 1.7. 
At z ~ 4 and 3, our results are in good agreement 
with the semi-analytic models of galaxy formation 
that have been added to the Millennium Simula- 
tion, but generally do not agree with the predic- 
tions at lower redshifts. 

3. The relatively small field of view of the KDFs im- 
poses a systematic dependence on the correlation 
lengths we infer due to the large integral constraint 
corrections that need to be applied. Howerver, by 
using artificial integral constraint corrections, we 
have shown that this dependence does not affect 
any of the important general trends that we ob- 
serve. 
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4. We find that the clustering strength of galaxies in 
all redshift samples is related to UV color. We pro- 
pose that this is due to the presence of the 2175 A 
absorption bump in more massive halos, which con- 
tain the older stellar populations and dust needed 
to produce the feature. This UV color relation pro- 
vides a caution against comparisons between sur- 
veys that have different selection criteria, and also 
for comparisons between different redshift samples. 

5. Our inclusion of fainter UV magnitudes has re- 
vealed additional similarities between Lyman break 
galaxies and -Bz-FT-selected galaxies at z ~ 2, such 
as 7?.-faint objects clustering more strongly than 1Z- 
bright ones, and correlation lengths greater than 10 
h^ 1 Mpc for subsamples of galaxies selected with 
each technique. 

6. Using the Millennium Simulation results, we de- 
rive halo occupation numbers for our samples and 
find that at z ~ 4 and 3, the results are generally 
consistent with approximately one galaxy per halo. 

7. At z ~ 2.2, we estimate the relation between stellar 
and dark matter mass, and compare it to model 
predictions. We find generally poor agreement, and 
believe this may be due to an underestimate of the 
stellar masses, whose estimates are dominated by 
ongoing star formation and likely miss the presence 
of older stellar populations. This is indeed what 
one would expect from galaxies in the process of 
turning down their star formation rates. 



In conclusion, we note that the addition of faint Ly- 
man break galaxies to clustering measurements suggests 
several interesting new trends; it will be important to 
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extend these findings with new surveys that reach the 
same faint magnitudes, but cover a larger area, in order 
to minimize the significance of systematics, in particular 
those due to the integral constraint. An additional future 
direction will be the color-selection of subsamples of the 
BM, BX, and LBG populations to investigate the effects 
of the shut-down of star formation in massive haloes. 
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